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Abstract — In  this  paper,  we  consider  the  problem  of  design¬ 
ing  transceiver  beamforming  vectors  for  multi-source  multi¬ 
destination  (MSMD)  wireless  networks  such  that  the  transmission 
power  of  each  source  is  minimized  while  the  signal-to-interference 
plus  noise  ratio  (SINR)  requirements  of  all  source-destination 
pairs  are  satisfied.  We  propose  an  efficient  iterative  algorithm 
to  design  the  transceiver  beamforming  vectors  and  address 
the  convergence  of  the  algorithm.  We  determine  a  necessary 
condition  as  well  as  a  sufficient  condition  for  the  algorithm 
to  converge  to  a  generalized  Nash  equilibrium  solution.  Espe¬ 
cially,  if  each  destination  has  only  one  antenna,  we  obtain  a 
necessary  and  sufficient  condition  for  the  algorithm  to  converge 
to  a  unique  generalized  Nash  equilibrium  solution.  Simulation 
results  show  that  the  proposed  iterative  algorithm  has  higher 
probability  of  convergence  compared  to  an  iterative  waterfilling 
(IWF)  approach.  For  example,  for  a  system  with  three  source- 
destination  pairs  and  SINR  requirement  of  2dB,  the  probability 
of  convergence  is  84%  with  the  proposed  algorithm  while  it  is 
only  66%  with  the  IWF  algorithm. 

Index  Terms:1  Multi-source  multi-destination  (MSMD)  wire¬ 
less  network,  MIMO  interference  channel,  transceiver  beamform¬ 
ing,  game  theory,  generalized  Nash  equilibrium. 

I.  Introduction 

Cognitive  radio  network  is  a  promising  technology  in  aim 
to  improve  spectrum  utilization  in  wireless  networks  where 
secondary  users  may  share  primary  users’  spectrum  resource 
without  causing  significant  interference  [1],  [2].  Cognitive 
radio  channel  model  can  be  characterized  as  interference 
channels  with  multi-sources  and  multi-destinations  (MSMD). 
To  improve  performance  of  interference  channels,  the  use 
of  multiple-input  multiple-output  (MIMO)  technology  turns 
out  to  be  an  effective  approach,  in  which  each  source  and 
destination  is  equipped  with  multiple  antennas.  In  recent  years, 
the  design  of  optimal  transceiver  beamformer  for  MIMO 
interference  channel  has  attracted  great  interest  [3]-[5],  where 
beamforming  vectors  were  designed  to  either  minimize  the 
total  power  consumption  of  the  system  or  maximize  Quality 
of  Service  (QoS)  of  the  system.  Various  QoS  criteria  have 
been  considered  in  the  literature  (see  for  example  [4],  [5], 
and  the  references  therein),  including  the  minimization  of  the 
weighted  sum  of  mean-square  error  (MSE),  the  maximization 
of  the  weighted  sum  of  user  data  rates,  and  the  minimization 
of  the  total  transmission  power  of  all  users. 

Game  theory  approach  [6]  has  been  utilized  to  design 
transceiver  beamformer  for  MSMD  networks  where  the  QoS 
of  each  source-destination  pair  is  maximized  subject  to  power 
constraint  of  each  source  [7]-[9].  In  [7],  an  iterative  algorithm 
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was  proposed  to  design  beamforming  vectors  and  it  was  shown 
that  the  algorithm  converges  to  a  Nash  Equilibrium  (NE)  point 
[6]  if  the  channel  matrices  are  of  full  column  rank.  In  [8],  a 
beamformer  design  algorithm  was  proposed  for  nonsingular 
channel  matrices  and  the  impact  of  multiuser  interference  on 
the  convergence  of  the  algorithm  was  discussed.  Furthermore, 
in  [9],  the  authors  extended  the  algorithm  to  accommodate 
arbitrary  channel  matrices  and  provided  a  unified  view  by 
proposing  an  iterative  waterfilling  (IWF)  algorithm. 

In  this  paper,  we  consider  the  problem  of  designing 
transceiver  beamforming  vectors  for  MSMD  MIMO  systems 
such  that  the  transmission  power  of  each  source  is  minimized 
while  the  signal-to-interference  plus  noise  ratio  (SINR)  re¬ 
quirement  of  each  source-destination  pair  is  satisfied.  We  know 
that  for  a  single-source  single-destination  MIMO  system,  the 
minimization  of  transmission  power  with  a  QoS  constraint  is 
equivalent  to  the  optimization  of  QoS  performance  with  a 
power  constraint  due  to  the  duality  of  the  two  optimization 
problems  [10].  However,  for  a  MSMD  MIMO  system,  the 
minimization  of  each  source-destination  pair’s  power  with  QoS 
constraints  is  no  longer  equivalent  to  the  optimization  of  QoS 
performance  with  power  constraints.  In  this  work,  we  propose 
an  iterative  algorithm  to  design  the  transceiver  beamforming 
vectors  which  has  a  higher  probability  of  convergence  than 
an  IWF  design  approach.  We  also  determine  a  necessary 
condition  as  well  as  a  sufficient  condition  for  the  algorithm 
to  converge  to  a  generalized  Nash  Equilibrium  solution.  Es¬ 
pecially,  if  each  destination  has  only  one  antenna,  a  necessary 
and  sufficient  condition  is  determined  for  the  algorithm  to 
converge  to  the  unique  generalized  Nash  Equilibrium  solution. 
Numerical  and  simulation  results  are  provided  to  illustrate  the 
proposed  algorithm  and  the  theoretical  development. 

Notation:  Uppercase  and  lowercase  bold  letters  denote 
matrices  and  vectors,  respectively.  (-)r,  (• )H  and  (•j1*  represent 
the  transpose,  the  Hermitian  and  the  Moore-Penrose  pseudoin¬ 
verse  of  a  matrix,  respectively,  p(-)  is  defined  as  the  spectral 
radius  of  a  matrix  [13].  A  matrix  B  or  a  vector  b  is  called 
non-negative  (positive)  if  each  component  of  the  matrix  or 
vector  is  non-negative  (positive),  denoted  as  B  >  0(>  0)  or 
b  >  0(>  0),  respectively. 

II.  System  Model  and  Problem  Formulation 

We  consider  a  MSMD  wireless  network  which  consists 
of  K  source-destination  pairs,  where  each  source  ,S). ,  k  = 
1,  2,  •  •  •  ,  K,  transmits  information  to  its  intended  destination 
Dk,k  =  1,2,  •••  ,K,  respectively.  We  assume  that  each 
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source  has  M  transmit  antennas  and  each  destination  has 
N  receive  antennas.  Let  be  the  transmitted  information 
symbols  at  source  Sk ■  Then  the  signal  sent  by  the  source  Sk 
can  be  expressed  as 

Sfc  =  w  kbk, 

where  Wfc  £  CMxl  is  a  beamforming  vector  at  source  Sk,  and 
D{|&fc|2}  =  1.  So,  the  average  transmission  power  at  source 

Sk  is 


£{||WfcMI  }  =  w£wfc.  (1) 

The  received  signal  at  destination  Dk  can  be  modeled  as 


K 


y  k  = 


Y. 

1=1, tyk 


nfc, 


(2) 


where  Hfcfc  £  C,VxM  is  the  channel  matrix  between  source 
Sk  and  destination  Dk,  H ik,  ( l  y  k),  is  the  channel  matrix 
between  source  Si  and  destination  Dk.  The  noise  vector  n,,; 
at  destination  Dk  is  assumed  to  be  an  A'-dimensional  white 
Gaussian  noise  vector  with  mean  zero  and  variance  crpjv- 
Without  loss  of  generality,  we  assume  ak  =  1,  1  <  k  <  K. 

At  each  destination  Dk,  a  receiver  beamforming  vector  u,  £ 
CNxl  is  applied  to  the  received  signal  yk,  then  the  resulting 
SINR  at  destination  Dk  can  be  represented  by 


wfHffcufeu^Hfcfcwfe 


SINR,  = 


ufufc  +  Ei=i ,i?k  wf  H^ufeuf  H;fcWi 


HttH, 


(3) 


in  which  the  interference  from  other  sources  is  treated  as 
additive  color  noise. 

For  any  given  SINR  requirement  (e.g.  QoS  constraint)  of 
each  source-destination  pair,  we  try  to  minimize  the  trans¬ 
mission  power  for  each  source.  The  optimization  problem  is 
formulated  as  follows: 


minimize  wjfwfc 
subject  to  SINRfc  >  7^ 


k  =  1,2,  ■  •  •  ,  K,  (4) 


where  7 k  is  the  SINR  requirement  of  the  source-destination 
pair  (Sk,  Dk).  Without  loss  of  generality,  we  assume  that  7 k  = 
1  for  all  k,  otherwise  we  can  always  obtain  an  equivalent 
problem  by  replacing  H kk  with  H kk/y/Tk-  According  to  (3), 
the  optimization  problem  can  be  written  as 


minimize 
subject  to 

„Ht 


wfwfe 


k  =  1, 2,  ■ 


.  K 


w"  GfcfcWfc  >  ufufe 


'  2si=i,ijtk  w; 


(5) 


H 


GikWi, 


A 


where  G ik  =  H/jJiifcu"  H,fc  for  1  <  l,  k  <  K. 


III.  Beamforming  Design  and  Convergence  of  the 
Algorithm 

A.  Transceiver  Beamformer  Design 

For  simplicity  of  notation,  let  us  denote  the  transmitter 
beamformer  set  as  W  =  (wi,...,Wa)  and  the  receiver 
beamformer  set  as  U  =  ( u  1 , ....  ua).  First,  we  consider 
receiver  beamforming  design  at  each  destination  for  any  fixed 


transmitter  beamformer  set  W.  Note  that,  at  each  destination 
Dk,  the  covariance  matrix  of  the  noise  plus  interference  is 
given  by 

K 

Rfc  =  I  n+  Y  (6) 

1=1, l^k 

which  can  be  estimated  at  destination  D,c  based  on  the  received 
signals.  The  receiver  beamforming  vector  u/r  that  maximizes 
SINRj.  at  destination  Dk  is  given  by  [4] 

ufc  =  R^HkkWk-  (7) 

The  resulting  maximal  SINK ,,  is 

SINRfc  =  wf  H^R-'HfcfcWfc.  (8) 


Now  let  us  consider  beamforming  design  at  the  transmitter 
side.  At  each  source  Sk,  with  given  receiver  beamforming 
vector  Ufc  and  fixed  other  users’  transmit  beamformer  vectors, 
to  maximize  Gkk^k  in  (5),  the  optimal  transmit  beam- 
former  w k  at  source  Sk  should  be  parallel  to  the  vector  II^.u,,.. 
Let  us  denote  the  normalized  vector  of  H^.u^  as  v,.  and 

A k  =  ujfHfcfcH^.Ufc,  then  we  have 

Dc  =  -iHffcUfe.  (9) 

V  Afc 

To  minimize  the  transmission  power  w;,/  vvy.  in  (5),  by  La- 
grangian  method,  the  optimal  transmit  beamformer  w k  should 
follow  with  the  equality  of  the  constraint,  i.e.,  wkGkk.y/k  = 
ufufc  +  Eti,^fe  wj1  Gii,  \\'i,  and  with  the  notation  in  (9),  the 
corresponding  transmission  power  w;,/  w,.  is  given  by 


1 


K 

IWfeir  =  ^(ufufc+  Y  wfG;fcw;). 

k  1=1, tyk 


(10) 


Therefore,  for  any  source  Sk,  1  <  k  <  K,  the  optimal  transmit 
beamformer  w,c  that  minimizes  its  transmission  power  subject 
to  the  QoS  constraints  should  satisfy  the  following  equation: 


wfc  = 


1/2 


G,fcW,) 


Vfc,  1  <  k  <  K. 


(ID 

According  to  (6),  and  (9),  the  optimal  transmit  beamformer 
w k  at  each  source  Sk  can  be  represented  as 


Wfc 


1  <  k  <  K, 


(12) 


in  which  the  covariance  matrix  Rfc  is  defined  in  (6)  that  can 
be  estimated  at  destination  Dk. 

From  (10)  and  (1 1),  we  can  see  that  the  power  minimization 
problem  can  be  formulated  as  a  fixed-point  problem  as  follows. 

Denote  x  =  (xi,X2,---  ,xk)T,  where  Xk  =  1 1 Wfc 1 1 2  for  any 
k  =  1,2,- ■■  ,  K ,  then  the  equation  (10)  can  be  written  as 

xk  =  ArOfcV-  +  EzIwfcvfGifcViz/)  for  any  1  <  k  < 
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K.  Let  us  define  a  mapping  F  :  R+  — >  R+  as,  F(x)  = 
(/i  (x) ,  h (x) ,  •  •  •  ,fK(x))T,  where 

1  K 

fk(x)  =  —  (ufufe  +  ^  rfGikYiXi),  1  <  k  <  K.  (13) 

fc  l-l.l-fl; 

Then,  to  find  the  optimal  transmit  beamformer  vectors 
(wi,  •  •  •  ,  w k)  that  satisfy  the  K  equations  in  (10)  is  equiva¬ 
lent  to  solve  the  following  fixed-point  problem: 


x  =  F(x). 


(14) 


A 


Furthermore,  let  us  denote  vector  a  =  (ai,  02, ax) 


where  ak  =  \k ,  and  matrix  B  with  components 

vfGifcV; 


R  ^ 

Okl  — 


^ k 

0, 


1  <  k  ^  l  <  A ; 
1  <  k  =  l  <  K, 


(15) 


then,  the  fixed-point  problem  in  (14)  can  be  written  as 


resulting  beamformer  set  (W,  U)  satisfy  (7)  and  (12)  for  each 
source-destination  pair,  i.e., 


Ufc 

Wfc 


Rfc  1HfcfcWfc, 
v/ uf  RfcUfe 


1  <  k  <  K, 

Hffcufc,  1  <  k  <  K. 


(21) 

(22) 


Note  that,  from  (21)  and  (22),  we  have 


ufc  RfcUfc 


uf  HfcfcWfc 

H 11 

U fc  H kk 


u"HfcfcHfefeufc 


HfcfcUfc 


uf  RfcUfc, 


which  implies  that  ujf  RfcUfc  =  1.  So,  the  iterative  process  in 
(22)  can  be  replaced  by 


Wfc 


1  <  k  <  K. 


(23) 


x  =  F(x)  =  a  +  Bx. 


(16) 


Note  that,  according  to  (9),  the  components  of  the  vector  a 
can  be  specified  as 


ak 


ufufc 

ufHfcfcHffeUfc  ’ 


1  <  k  <  K, 


(17) 


and  the  components  of  the  matrix  B  can  be  specified  as 


ufHnH/fuzufHfcfcHffcUfc’ 


1  <k^l<K.  (18) 


We  observe  that  if  I  B  is  invertible,  the  fixed-point  problem 
in  (16)  has  a  unique  solution  which  is  given  by2 


x=  (I  —  B)_1a. 


(19) 


If  the  solution  x  is  positive,  then  the  optimal  transmitter 
beamforming  vectors  (wi,  •  •  •  .  w^-j  are  determined  as 

Wfc  =  V^fcVfc,  1  <  k  <  K,  (20) 

where  Vfc  is  specified  in  (9),  and  the  corresponding  optimal 
receiver  beamforming  vectors  (ui,  •  •  •  ,  u  k)  can  be  determined 
based  on  (7)  accordingly. 

We  may  view  the  MSMD  beamformer  design  problem  as 
a  noncooperative  Nash  equilibrium  game  [6],  in  which  each 
source-destination  pair  is  an  active  player  and  they  compete 
each  other  to  minimize  their  own  transmission  power.  The 
payoff  of  each  player  is  its  own  transmission  power  and  the 
strategy  of  each  player  is  to  act  aggressively  and  selfishly  to 
adjust  its  transceiver  beamformer  vectors  Ufc  and  Wfc  according 
to  (7)  and  (12),  respectively.  We  can  see  that  the  strategy  of 
each  player  depends  on  others’,  so  it  is  a  generalized  Nash 
equilibrium  (GNE)  game  [6]  and  the  corresponding  equilib¬ 
rium  is  called  a  GNE  solution.  A  transceiver  beamformer  set 
(W,  U)  is  a  GNE  solution  if  for  each  source-destination  pair 
(S'fcjDfc),  given  the  other  pairs’  transceiver  beamformers,  the 
transmission  power  of  source  Sk  cannot  be  further  reduced  by 
optimizing  the  transceiver  beamformer  design  of  this  pair.  The 


2Note  that  if  I  B  is  not  invertible,  we  may  use  the  Moore-Penrose 
pseudoinverse  (I  —  B)**  to  represent  the  solution  x  in  (19). 


In  the  following,  we  propose  an  iterative  algorithm  to 
design  transceiver  beamformer  set  (W,  U)  that  minimizes 
each  source’s  transmission  power  while  satisfying  the  SINR 
requirements  of  all  source-destination  pairs. 

Algorithm: 

•  Step  1 :  Initialize  the  transceiver  beamformer  set  (W,  U) 

with  random  vectors; 

•  Step  2:  Determine  matrix  B  according  to  (18),  and  set 

x  =  (I  —  B)**a,  where  a  is  given  by  (17); 

•  Step  3:  if  x  >  0,  then 

update  wj”"1  using  (20)  for  all  1  <  k  <  K; 
else 

update  w[r^  using  (12)  for  k  with  Xk  <  0; 
end; 

.  Step  4:  if  ||w^"+1)  —  w^H2  <  e  for  all  1  <  k  <  K, 
then 

output  the  transceiver  beamformer  vectors 
and  for  all  1  <  k  <  K. 

else 

update  uf'l>  by  (7),  and  goto  Step  2; 
end; 

•  Step  5:  if  the  maximal  number  of  iterations  reaches, 

then 

output  failure  message  and  stop; 
else  goto  Step  1. 

The  convergence  of  the  above  iterative  algorithm  is  ad¬ 
dressed  in  the  next  subsection.  We  note  that  an  alternative 
iterative  waterfilling  (IWF)  approach  [9]  may  be  utilized  to 
design  the  transmit  beamformer  set  W.  In  this  case,  the 
transmit  beamformer  for  each  source  Wfc  is  chosen  as  the 
eigenvector  of  Hfcfc  corresponding  to  the  maximal 

eigenvalue  \max,  which  results  in  w^Wfc  =  l/\max-  As  we 
will  see  in  the  simulation  results,  our  proposed  algorithm  has 
a  higher  probability  of  converging  to  a  GNE  solution  than  the 
IWF  approach. 
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B.  Condition  to  the  Convergence  of  the  Algorithm 

In  this  subsection,  we  determine  a  necessary  condition  as 
well  as  a  sufficient  condition  for  the  convergence  of  the 
proposed  iterative  algorithm.  We  first  develop  two  lemmas 
which  are  critical  to  obtain  the  main  result. 

Lemma  1:  For  a  matrix  B  >  0  and  a  vector  a  >  0,  the 
necessary  and  sufficient  condition  for  equation  (I  —  B)x  =  a 
to  have  a  unique  and  non-negative  solution  is  that  the  spectral 
radius  p(B)  is  less  than  1,  where  p(B)  =  max{\Ai\}  and  A,’s 
are  eigenvalues  of  B. 

Proof:  If  p(B)  <  1,  we  know  that  both  B  and  I  —  B  are 
invertible  [13].  Moreover,  since  (I  — B)'1  -  ^(J  Bn,  which 

is  non-negative  with  B  >  0,  thus  with  positive  a,  the  equation 
(I  —  B)x  =  a  has  a  solution  x  =  (I  —  B)-1a  which  is  unique 
and  non-negative. 

On  the  other  hand,  for  any  non-negative  matrix  B,  according 
to  the  Perron-Frobenis  theorem  [13],  there  exists  a  non¬ 
negative  vector  d  /  0  such  that  BTd  =  p(B)d.  So,  if  the 
equation  (I  —  B)x  =  a  has  a  non-negative  solution,  we  have 

(1  -  p(B))dTx  =  dT(I  -  B)x  =  dTa  >  0. 

Since  dTx  >  0,  so  1  —  p(B)  must  be  positive.  Thus,  we  have 
the  necessary  condition  p(B)  <1.  ■ 

Lemma  2:  If  a  matrix  B  £  RAxA  has  the  following  form: 

R  _  f  bu  >  0,  1  <  k  ^  l  <  K; 

~  X  bkk  =  0,  1  <  k  <  K, 

then  the  spectral  radius  of  B  is  less  than  1,  i.e,  (AW)  <  1,  if 
and  only  if  there  exist  K  positive  values  X2,  •  •  -  ,  x k  such 
that 

1  K 

—  Y  WklXl  <1,  1  <  k  <  K.  (24) 

1= l,ljik 

Proof:  If  (AW)  <  1,  then  by  Lemma  1,  for  any  vector 
a  >  0,  there  exists  a  unique  and  non-negative  vector  x  = 
(xi,X2,  •  •  •  , Xk)  that  satisfies  the  equation  (I  —  B)x  =  a.  In 
this  case,  we  have  (I  —  B)x  >  0,  which  leads  to  the  condition 
in  (24).  Conversely,  if  the  condition  in  (24)  holds,  which  is 
equivalent  to  (I  —  B)x  >  0,  then  there  exists  a  positive  vector 
a  such  that  (I  —  B)x  =  a,  which  means  that  the  equation  has 
a  non-negative  solution  x.  According  to  Lemma  1  again,  we 
have  (AW)  <  1,  which  completes  the  proof.  ■ 

Based  on  Lemmas  1  and  2,  we  are  able  to  determine  in  the 
following  a  necessary  condition  for  the  proposed  algorithm  to 
converge  to  a  GNE  solution. 

Theorem  1:  A  necessary  condition  for  the  proposed  algo¬ 
rithm  to  converge  to  a  GNE  solution  is  that  there  exists  a 
vector  set  {u^.  £  CNxl  :  ||ufc||2  =  1, 1  <  k  <  K}  such  that 

p(CD2)  <  1,  (25) 


where  the  matrices  C  and  D  are  given  by 


r  - 


D  kl  = 


U k 

0, 

i 

0, 


i  1  <  k  l  <  K) 
1  <  k  =  l  <  K, 

1  <  k  =  l  <  K\ 
l<k^l<K. 


(26) 

(27) 


Proof:  If  a  GNE  solution  exists,  it  means  that  there  exists 
a  transceiver  beamforming  vector  set  (W,  U)  satisfying  the  two 
recursive  equations  (21)  and  (22).  Furthermore,  it  implies  that 
the  fixed-point  problem  in  (16)  has  a  solution  of  xk  =  \\[!  wk , 
1  <  k  <  K,  where  vvy.  is  given  in  (22)  or  equivalently  in 
(23).  Following  the  expression  in  (23),  the  solution  xk  can  be 
specified  as 


Xk 


l<k<  K. 


(28) 


By  substituting  the  above  solution  into  the  equation  in  (16), 
we  have 


h  ^ufHnH^ufeufH;feHAu, 

ufc  Ufc+;  - 


(u;  H/;H//u;) 


1,  1  <  k  <  K.  (29) 


Denote  q  =  (q1, ...,  qK)T ,  where  qk  =  }  and  let  ufc  = 

uk 

^ /qf-Uk  for  1  <  k  <  K,  then  ||fife ||2  =  1,  and  with  the  matrices 
defined  in  (26)  and  (27),  the  equation  in  (29)  can  be  written 
as 

1+  Y  ^nSL  =  i,  i  <k<K,  (30) 

qk  qk 

which  leads  to  a  matrix  representation  as  follows: 

(I  —  CD2)q  =  1.  (31) 


According  to  Lemma  1,  the  existence  of  the  above  equation 
implies  that  the  spectrum  radius  of  the  matrix  CD2  must  be 
less  than  1,  i.e.  p(CD2)  <  1.  Thus,  we  prove  the  theorem 
completely.  ■ 

For  some  channel  conditions,  if  the  necessary  condition  in 
Theorem  1  is  not  satisfied,  the  proposed  iterative  algorithm 
will  diverge.  For  example,  let  us  consider  a  system  with  two 
source-destination  pairs  ( K  =  2)  with  M  —  2  and  N  =  1.  If 
the  channels  are  given  by 


Hn 

H2i 

Hi2 

h22 


(0.1078  +  0.7184*, 
(-3.1297-  1.6945i, 
(—4.2857+  1.8133i, 
(1.0186+  1.0264*, 


0.1189 +  0.739H), 
-1.2949  +  3.0532*), 
4.0810  -  3.0028*), 
0.2113  +  0.1098*), 


then  for  any  unit  vector  u k,k  =  1,2,  the  matrix  CD2  is 


CD2 


0  21.7475  \ 

22.0352  0  J  ' 


(32) 


The  spectral  radius  of  the  matrix  CD2  is  p(CD2)  =  21.8909 
which  is  larger  than  1.  In  this  case,  both  the  proposed 
algorithm  and  the  IWF  algorithm  diverge  as  shown  in  Fig.l. 

In  the  following,  we  obtain  a  sufficient  condition  for  the 
proposed  algorithm  to  converge  to  a  GNE  solution.  Before  we 
present  the  theorem,  we  introduce  a  lemma  and  its  proof  is 
omitted  due  to  space  limitation. 

Lemma  3:  For  a  K  x  K  matrix  B  >  0  and  a  K  di¬ 
mension  vector  a  >  0,  if  there  is  a  non-negative  vector 
x  =  (xi,  X2, Xk )T  such  that  (I  —  B)x  >  a,  then  there 

exists  a  non-negative  vector  x  =  (xi,X2 ,---,Xk)T  such  that 
(I  —  B)x  =  a  and  xk  <  xk  for  all  1  <  k  <  K. 
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Theorem  2:  The  proposed  algorithm  converges  to  a  GNE 
solution,  if  for  any  vector  set  {u*  G  CNxl  :  ||ua;||2  =  1,1  < 
k  <  K},  the  following  condition  is  satisfied 

p(CD2)  <  1,  (33) 


where  the  matrices  C  and  D  are  defined  in  (26)  and  (27), 
respectively. 

Proof:  For  any  vector  set  {u*  £  CNxl  :  ||ufc||2  = 
1,1  <  k  <  K } ,  we  observe  that  the  matrix  B  in  (18)  can 
be  represented  as  B  =  DCD,  where  the  matrices  C  and  D  are 
defined  as  in  (26)  and  (27),  respectively.  So,  we  have 


P(B)  =  p(DCD)  =  p(CD2).  (34) 


If  p(CD2)  <  1,  then  p(B)  <  1.  According  to  Lemma 
1,  the  fixed-point  problem  in  (16)  always  guarantees  a 
non-negative  solution  for  any  receiver  beamforming  vec¬ 
tor  set  U.  Let  us  denote  =  (u$"\  u^, ...,  u^)  and 

W(n)  =  (wi”\w2n\...,w^)>  with  uaT^  anc*  wLn)  denot¬ 
ing  the  transceiver  beamforming  vectors  of  the  /c-th  source- 
destination  pair  at  n  iteration  of  the  proposed  algorithm. 
At  any  iteration  n,  the  transceiver  beamforming  vector  set 
(W(n),U(n))  is  chosen  to  satisfy  the  equations  in  (16),  or 
equivalently  in  (10).  It  implies,  from  (9)  and  (20),  that 
the  SINR  of  the  k- th  source-destination  pair  given  by  the 
transceiver  beamforming  vector  set  is  equal  to 

its  minimal  requirement,  i.e.  SINR*(W^,U^)  =  1  for  all 

I  <  k  <  K.  Atn  +  1  iteration,  receiver  beamformer 

is  chosen  to  maximize  SINR*  for  given  transmit  beamformer 
W(n),  so  SINRfc(W(n),U(n  +1i)  >  1,  which  means  there  exists 

a  transmit  beamformer  set  W  =  (\v'  j .  #9 . W/<-)  such  that 

SINRfc(W,U(n+1))  =  1  for  all  1  <  k  <  K  and  ||wfc||2  < 

I I  w*n)  1 12 ,  for  all  1  <  k  <  K.  For  given  receiver  beamformer 

set  U*-"-1"1),  since  is  the  optimal  transmit  beamformer 

set  which  minimizes  the  transmission  power  ||wfc||2  for  all 
1  <  k  <  K,  so  ||w^”+1^||2  <  1 1 Wfc 1 1 2  for  all  1  <  k  <  K. 
Therefore  we  have  ||wj,n+1)||2  <  ||w["')||2  for  all  1  <  k  <  K. 
It  means  that  the  transmission  power  ||wfc||2,  which  has  lower 
bound  0,  decreases  at  every  iterative  step  for  all  1  <  k  <  K, 
thus  the  proposed  algorithm  converges  to  a  GNE  solution.  ■ 


Interestingly,  based  on  Theorems  1  and  2,  we  are  able  to 
obtain  a  necessary  and  sufficient  condition  for  the  convergence 
of  the  iterative  algorithm  if  each  destination  has  only  one 
antenna  ( N  =  1).  In  this  case,  each  receiver  beamformer  is  a 
scalar,  i.e.  u*  =  1  for  1  <  k  <  K,  and  the  matrices  in  (26) 
and  (27)  are  reduced  to: 


r  \\HlkHjj 
kl  0, 


Da-z  = 


0, 


1  <  k  y  l  <  K] 
1  <k  =  l<K, 

1  <k  =  l<I<- 
1  <  k  y  l  <  K , 


(35) 

(36) 


which  depend  only  on  the  channel  conditions.  Therefore, 
according  to  Theorems  1  and  2,  the  proposed  iterative  al¬ 
gorithm  converges  to  a  unique  GNE  solution  if  and  only  if 


p(CD2)  <  1.  We  summarize  the  above  discussion  in  the 
following  result: 

Corollary:  If  each  destination  has  only  one  antenna  (N  = 
1),  then  the  proposed  iterative  algorithm  converges  to  a  unique 
GNE  solution  if  and  only  if 

p(CD2)  <  1, 

where  the  matrices  C  and  D  are  specified  in  (35)  and  (36), 
respectively. 

IV.  Numerical  and  Simulation  Results 

In  this  section,  we  carry  out  some  numerical  and  simulation 
studies  to  illustrate  the  performance  of  the  proposed  iterative 
algorithm  and  compare  it  with  the  IWF  algorithm.  We  assume 
that  there  are  three  source-destination  pairs  ( K  =  3),  each 
source  has  two  transmit  antennas  (M  =  2)  and  each  destina¬ 
tion  has  two  receive  antennas  (N  =  2).  Channel  coefficients 
are  assumed  to  be  Rayleigh  distributed  with  unit  variance.  The 
SINR  requirement  of  each  source-destination  pair  is  assumed 
to  be  7*  =  1  (1  <  k  <  K). 

In  Figs.  2  and  3,  we  illustrate  the  performance  of  the  pro¬ 
posed  iterative  algorithm  with  a  randomly  generated  channel 
condition.  Specifically,  Fig.  2  shows  the  required  transmission 
power  at  each  source  in  each  iteration.  We  observe  that  with 
the  proposed  algorithm,  the  transmission  power  of  each  source 
converges  after  10  iterations,  which  leads  to  a  GNE  solution. 
With  the  IWF  algorithm,  however,  the  transmission  power  of 
each  source  does  not  converge  and  fluctuates  periodically. 
Fig.  3  plots  the  resulting  SINR  of  each  source-destination 
pair.  We  can  see  that  with  the  proposed  algorithm,  the  SINR 
requirement  of  each  source-destination  pair  is  satisfied  after  10 
iterations.  But  with  the  IWF  algorithm,  the  SINR  requirement 
of  the  three  pairs  cannot  be  satisfied  at  the  same  time  and  the 
resulting  SINRs  fluctuate  in  different  iterations. 

In  Fig.  4,  we  compare  the  probability  of  convergence  for 
the  proposed  algorithm  and  the  IWF  algorithm  with  different 
SINR  requirements  71  =  72  =  73  £  (OdB,  15dB)  over  2000 
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Fig.  2.  Transmission  power  at  each  source  with  the  proposed  algorithm  and 
the  IWF  algorithm,  K  =  3,  M  =  N  =  2. 


Fig.  3.  Resulting  SINR  at  each  destination  with  the  proposed  algorithm  and 
the  IWF  algorithm,  K  =  3,  M  =  N  =  2. 


randomly  generated  channel  realizations.  We  can  see  that  the 
convergence  probability  of  the  proposed  algorithm  is  much 
higher  than  that  of  the  IWF  algorithm.  For  example,  with  a 
SINR  requirement  of  2dB,  the  probability  of  convergence  is 
66%  by  the  IWF  algorithm  while  it  is  84%  by  the  proposed 
algorithm.  With  a  SINR  requirement  of  5dB,  the  probability  of 
convergence  is  improved  from  36%  to  66%  with  the  proposed 
algorithm. 

V.  Conclusion 

In  this  paper,  we  designed  transceiver  beamforming  vec¬ 
tors  for  MSMD  MIMO  wireless  networks  such  that  the 
transmission  power  of  each  source  is  minimized  while  the 
SINR  requirements  of  all  source-destination  pairs  are  satisfied. 
We  proposed  an  iterative  algorithm  to  design  the  transceiver 
beamforming  vectors  and  discussed  the  convergence  of  the 
algorithm  which  depends  on  channel  condition.  We  determined 
a  necessary  condition  as  well  as  a  sufficient  condition  for  the 
algorithm  to  converge  to  a  GNE  solution.  If  each  destination 


Fig.  4.  Probability  of  convergence  of  the  proposed  algorithm  and  the  IWF 
algorithm,  K  =  3,  M  =  N  =  2. 


has  only  one  antenna,  we  obtained  a  necessary  and  sufficient 
condition  for  the  algorithm  to  converge  to  a  GNE  solution. 
Simulation  results  show  that  the  proposed  iterative  algorithm 
has  higher  probability  of  convergence  compared  to  the  IWF 
algorithm.  For  example,  for  a  system  with  three  source- 
destination  pairs  and  SINR  requirement  of  2dB,  the  probability 
of  convergence  is  84%  by  the  proposed  algorithm  while  it  is 
only  66%  by  the  IWF  algorithm 
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